Upper critical field in a spin-charge separated superconductor. 
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It is demonstrated that the spatial decay of the pair propagator in a Luttinger liquid with spin charge 
separation contains a logarithmic correction relative to the free fermi gas result in a finite interval 
between the spin and charge thermal lengths. It is argued that similar effects can be expected in 
higher dimensional systems with spin charge separation and that the temperature dependence of 
the upper critical field Hc2 curve is a probe of this effect. 
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Recent experiments^'^ on low transition temperature 
single layer cuprates indicate that the upper critical field 
curve in these materials has an anomalous upward cur- 
vature extending to temperatures below a small fraction 
of Tc- This behavior contrasts strongly with the mean 
field or weak coupling BCS result which gives an approx- 
imately parabolic shape for the Hc2 curve. This discrep- 
ancy is especially interesting because weak coupling the- 
ory could be expected to apply to systems with low Tc-^ 

Near to Tc{H) the Hnearized gap equation becomes an 
homogeneous integral equation with a kernel which rep- 
resents the normal state electron pair propagator in a 
magnetic field. Gorkov'' included the effect of the gauge 
dependent phase acquired by a propagating electron, but 
assumed that the role of Landau quantization^ could be 
neglected due to temperature or lifetime effects. Special- 
izing to two dimensions, the linearized Gorkov equation 
in the symmetric gauge becomes: 

A(r) =g j dPr'Koir', /3) eM^-^)Hr + r') (1) 

where A(r) is the space dependent energy gap, ji is the 
inverse temperature, g is the (assumed field independent) 
pairing interaction. The magnetic length I is related to 
the appHed field by -H" = 4io{2itI'^)~^ where <ji»o is the fiux 
quantum hc/e. Kci{r) is the fermion pair propagator or 
pair susceptibility in real space, in the absence of the 
external field and the pairing interaction g. Rajagopal 
and Vasudevan® pointed out the existence of a solution of 
Eq. 1 of the form A(r) = Aoe^^li-) . When substituted 
into Eq.l this leads to the pairing instability condition 
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where ro is a lower cutoff. Ko{l,l3) is the uniform pair 
susceptibility in the presence of an external field. 

Eq.2 is the mean-field instability condition defining the 
critical magnetic length lc{l3). It states that the upper 
critical field is a contour of constant pair susceptibility 



in the field-temperature plane. The temperature depen- 
dence of Hc2 depends on the details of the normal state 
pair propagator. In the case of free fermions in D dimen- 
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When r < fivf = ^, the real space pair susceptibility 

decays as a power law A'q ^ oc r^^ in 2D. At distances 
longer than the thermal length ^ the pair propagator is 
exponentially small. Thus the integral Eq. 2 depends 
logarithmically on a long distance cutoff which is equal 
to ^ or /, whichever is shorter. Substitution of Eq. 3 into 
Eq.2 leads to the conventional 2D almost paraboHc Hc2 



curve. 



We can explore the consequences of an anomalous 



pair susceptibility via the ansatz A'o(r, /?) 



'Z?-" for 



ro < r < I. For a 2D Eermi liquid we have seen that 
fj, = 2 and f = 0. The low temperature part of the Hc2 
curve (/ -C £,) is determined by this scaling. A weakly ro 
dependent result requires yU < 2 and, with this assump- 
tion, Eq.2 gives Ic « /J^-c or Hc2 ~ T^-c. A power 
law "divergence" of Hc2 at low temperature occurs when 
V <0. 

The exponents v,/! can be directly related to scal- 
ing properties of the normal state single particle Green 
function g{k,ui) provided it is assumed that the ver- 
tex function has a trivial scaling. That is, we assume 
(w„ = 27m//3), 



A'o(r,/3)=/3- 
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Generalising an argument due to Balatsky''', g{Ak, Aui) = 
A-"g{k,uj) implies that A'o(r,/3) = /3"''"^"/(f ). Thus 
our scaling ansatz for Kq (r, (3) is recovered provided that 
f{x) ~ x^^ which leads to 
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FIG. 1. The (/-dependent pair susceptibility in a one di- 
mensional Luttinger model with §4 interaction only, whose 
Fourier transform is used in the text. Note that there are no 
vertex corrections in the absence of 172 coupling between left 
and right branches. 

Since // < 2, divergence of the upper critical field at low 
temperature is possible if a < —1.* 

While the experimental Hc2 curves show pronounced 
upward curvatures, the present data appear to be consis- 
tent with a finite upper critical field at zero temperature 
(i.e. V = 0). This suggests we pay special attention to 
the case of Fermi- liquid type scaling a = —1. In addi- 
tion to Fermi liquids, one-dimensional spinful Luttinger 
models with finite §4 interaction only satisfy a = —1. §4 
represents forward scattering of opposite spin electrons 
near the same fermi point. This produces spin charge 
separation but no anomalous scaling. For example, Z, 
the discontinuity in the fermion momentum distribution 
function at kf , equals unity. It is natural to ask how the 
spatial dependence of the pair propagator is changed in 
this case. Note that microscopic models such as the ID 
large U Hubbard model give small but finite values of §2 , 
the forward scattering interaction between left and right 
movers. 

The form of the ID spin-charge separated retarded 
Green function in real space and time for right and 
left branches is'^-iO; g^{r,t) = -^@{t){[0+ + i{x ^ 
v,t)]-^[G+ + i(x ^F Vct)]-i + (x ^ -x,t ^ -i)}e±^*/'' 
where the spin and charge velocities are Vc,s = Vfil^^). 
In momentum and frequency space, this is^° 
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The pair susceptibility (shown in Fig.l) can be expressed 



as. 



Ko{r,l3) = - dujtanh(Puj/2)Re{A{r,uj)B{-r,-uj)} 

(7) 

where A{r,uj) = -^ J dkexp{ikr)lm.g~^{k,uj) and 
B{r,ij) = -^ J dkex-p{ikr)Reg~^{k,uj). Note that there is 
no vertex function in Eq.7 in the absence of §2 coupling 
between left and right hand branches. We find that 
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and B{r,io) = isgn{r)A{r,uj). For the free fermion case 
{vs = Vc), A{r, to) has an oscillatory, non-decaying behav- 
ior. When Vs 7^ Vc the Bessel function leads to an addi- 
tional r^ 2 decay beyond a distance « VsVc{vc—Vs)^^u>^^ . 
This distance becomes large as the fermi surface is ap- 
proached. 

Eqs. 7 and 8 give 
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where we have introduced the effective thermal length 
^ = (3v where v~^ = v~^ + v~^ and the dimensionless 
parameter 
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The D = 1 free fermion result Eq.3, is recovered upon 
setting A = and v = Vf/2}^ 

The uniform pair susceptibility is obtained by integrat- 
ing over r with short distance cutoff ro -C f. Alterna- 
tively a high frequency cutoff can be introduced in Eq.9. 
Adopting the latter procedure, we find^^ 
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where K is the complete Elliptic function of the first 

kind. As A ^ i, K{2\) -^ -^ log(l - 2\) k \ log (f^) . 

Thus §4 interactions enhance the pair susceptibility rel- 
ative to the non- interacting result in this limit. 

We have not succeeded in evaluating the pair propaga- 
tor Eq.9 for general A in closed form. However progress 
can be made in the limit Vs -C Wc , where departures from 
Fermi liquid behavior are most significant. In this limit 
the dominant contribution can be obtained using the sub- 
stitution sin3;Jo(Aa;)^ -^ (7ra;)~^ cos(l — 2\)x which is 
valid under the integral sign in Eq.9. Thus, in terms of 
spin (^s = fivs) and charge (^c = (ivc) thermal lengths, 
we find 
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The essence of this result is the opening up of a new 
regime in the window ^s -C r -C ^c- In this regime 
the pair susceptibility decays more slowly than the free 
fermion result r~^ by a factor log(^c/^)- For r -C ^s,£,c, 
the behavior is r^^ with an enhancement factor log ( 7^ ) , 
while for r ^ (3vs,(3vc, A'o(r, /?) is exponentially small. 
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FIG. 2. Temperature dependence of the upper critical field 
on a normalized plot for (a) Vs = Vc, (b) v^ = O.lvc, (c) 
Vs = O.Ol-Uc, (d) Vs = O.OOl-Uc. 

We now ask whether an anomaly of this type, if present 
in a 2D system, would have an observable effect on Hc2- 
For example, if, following Anderson^'', we take the Green 
function Eq.6 as a model in D = 2, and continue to 
assume that vertex corrections can be neglected, then 
A'o(r, /3) has the form Eq.l2 with the prefactor r^^ re- 
placed by r^2 15 -y^g then solve Eq.2 for the critical field. 
The principle effects are: (1) The relationship between 
the zero temperature critical field and Tc is altered com- 
pared to the Fermi liquid result. When Vg = Vc, this is 



Vg -C Vc we find 
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(2) -H"c2(0) is enhanced relative to the slope of the critical 
field curve near to T^. When w, = w^ it is well known 



that Hc2{0) « Tc 
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.13 When Vg -C Vc we find 



^^logfe) 



dHr 



dT 



(14) 



Tc 



Numerical solution for dependence of Hc2{TY^ is shown 
in Fig. 2. As indicated above, spin-charge separation 
leads to upward curvature of the Hc2 line. Saturation 
behavior of the upper critical field always occurs at suf- 
ficiently low temperatures, when I -C ^c,£,s- Conversely, 
linear variation of Hc2 with temperature occurs near zero 

field />^,, 6- 

In summary, the upper critical field is a probe of both 
spatial and temperature (frequency) dependence of nor- 
mal state pair correlations. In a BCS superconductor 
the upper critical field satisfies the approximate relation- 
ship av'i{H I 4>q) + T"^ = Tc where a is a numerical con- 
stant. This result is a direct manifestation of Fermi liq- 
uid theory for the normal state, which is characterized 



by a single quasiparticle velocity scale Vf. The effect of 
spin charge separation is visible in the upper critical field 
through the presence of distinct thermal lengthscales for 
spin and charge. 
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